Abstract. Using a numerical simulation of compressible 
convection with radiative transfer mimicking the solar 
photosphere, we compare the velocity field derived from 
granule motions to the actual velocity field of the plasma. 
We thus test the idea that granules may be used to trace 
large-scale velocity fields at the sun's surface. Our results 
show that this is indeed the case provided the scale sepa- 
ration is sufficient. We thus estimate that neither velocity 
fields at scales less than 2500 km nor time evolution at 
' scales shorter than 0.5 hr can be faithfully described by 
. granules. At larger scales the granular motions correlate 
' linearly with the underlying fluid motions with a slope of 
i; < 2 reaching correlation coefficients up to ~ 0.9. 

^ , Key words: Convection - Sun: granulation - Sun: pho- 
' tosphere 
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1. Introduction 



Since the work of |November fc Simon (1988^ granules and radiative transfer^ originally developed by two 



have V^Ti ir. ir-cr-a hr.r-;..»rif al flr^wc at iha cn.r-for.^ of th e authors (for a description see [Stein fc Nordlund 



of the Km, r\Hmp]y mesnsralp flows Find siipprgraniilatinn 



([November, 1989|; jStrous, 19954 [fc|; [Roudier et al., 1999| ; 
Rieutord et al., 20001 ). However, assuming that granules 



behave like passive scalars is a rather strong assumption 
regarding the nature of the granules: these are dynam- 
ical structures which are far from being passive. Unfor- 
tunately, the validity of this assumption has never been 
assessed and one relies on the hope that advection of gran- 
ules (as intensity structures) is statistically dominant com- 
pared to noise processes like the diffusion of temperature 
fluctuations and small-scale motions induced by granules. 

To make further progress on this issue, we used a sim- 
ulation of convection at the sun's surface to test the track- 
ing properties of granules. The simulation provides an ob- 
servable - a time series of two dimensional images of the 
emergent intensity - together with the underlying three 
dimensional velocity and temperature fields. This allows 
to test and compare horizontal velocities as measured by 
different granule tracking techniques against the actual 
flow velocities. The tracking techniques are presently two, 
namely LCT, for Local Correlation Tracking, and CST 
for Coherent Structures Tracking. LCT has been devel- 



oped by November and collaborators (see November & 



Simon, 198S) and determines the flows from an optimiza- 



tion of the correlation between two subimages belonging 
to successive images. CST was first proposed by ^trous 



(1995a, E|) and recently dev eloped by Roudier et al. (1999[ ) 



and Rieutord et al. (2001); it decomposes each image of 



the solar surface into a set of granules whose trajectories 
are used to derive the velocity fields. Both methods have 
shortcomings whose effects can be quantified by the above 
comparison. 

The simulation used in this letter has been performed 
using a compressible convection code coupling fluid mo- 



1998). For this particular run, which aims at simulating 



Send offprint requests to: M. Rieutord 



supergranulation, certain trade-offs have been made be- 
tween physical realism and computational demands: the 
radiative transfer was treated in grey approximation (with 
frequency independent opacity, with a dependence on tem- 
perature and pressure similar to that of the solar contin- 
uum opacity) and the horizontal resolution was chosen to 
be a rather coarse 95 km. The restrictions made it afford- 
able to study a large volume (30 x 30 Mm^ wide and 3 Mm 
deep represented by 315x315x82 grid points) which con- 
tained several hundred granules at any given instant in 
time. 

In section 2 of this letter we shall present a global 
view of the velocity field, following Euler's viewpoint while 
in section 3 we try to characterize the granules in their 
ability at tracing the flow field, thus adopting Lagrange's 
viewpoint. Our conclusion is that granules are able to trace 
statistically the large-scale flows but lead to a systematic 
underestimation of the actual velocities. 

2. Euler's view 

As granules are extended test particles, the measured ve- 
locity fields are much less resolved than the one issued 
from the simulation. Typically, our granule tracking tech- 
nique yields a velocity field on a 45x45 grid, i.e. seven 
times coarser than the original 315x315 pixels. Hence, for 
comparison, simulated velocity fields are rebinned (aver- 
aged) to this coarser resolution. Granules also decrease 
the time resolution and velocity fields issued from granule 
tracking are usually averaged over a time window longer 
than 5 min; here, we shall consider three time-windows 
with durations of 1000 s, 1 h, and 2 h. 

But granules are also three-dimensional structures and 
therefore they "feel" the large-scale velocity fields aver- 
aged over some range of depths. The approximate depth 
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Fig. 1. Correlations between the velocity field obtained from granule tracking and the actual velocity field as a function 
of depth and time-window; both methods (LCT and CST) give very similar results. A layer 78 km thick has been 
used. 
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Fig. 4. For three length scales (2", 2^, 2^), we represent the actual velocity {v which is either Vx or Vy) as a function 
of the measured velocity u for a two-hours average. The clouds represent the 45^ grid points. The best fits (dashed 
hne) show that granules underestimate actual velocities by a factor which varies from 2.1 at small scales (2°) to 1.6 
at large scales (2^). 



and thickness of the contributing layer need to be deter- 
mined. 

Using the three above mentioned time-windows, we 
plotted in Fig. ^ the linear correlation between veloc- 
ity fields issued from granule tracking u (for a descrip- 
tion of how this field is derived, see Roudier et al. 1999 
or Rieutord et al. 2001) and the "original" ones v &s a. 
function of depth; this correlation is defined by Cy = 
{u ■ v) / y/ {u^) {v^). This figure clearly shows that the cor- 
relation increases with the length of the time window 
hence showing that granules are best for tracing large- 



scale flows which evolve on long time scales. It also shows 
that the correlation is best at a depth below the t — 1 
surface (here 2=0); this is clearly emphasized by Fig. || 
where we see that the optimal depth (where the correla- 
tion is maximum) is between 200 km and 300 km. 

We also tested the dependence of correlations with re- 
spect to the thickness of the layer and found that it is 
weak: variations of correlations are of 2 or 3% when the 
thickness of the layer is varied between 40 and 900 km. 

The foregoing results show that granules do trace long 
time-averaged flows, thus we should observe a better cor- 
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Fig. 2. Optimal depth for different thicknesses (line type) 
of the contributing layer as a function of the time- window. 
We see that the depth at which correlation is maximum, 
remains around 250±30 km. 



relation when small spatial scales are filtered out. This is 
indeed the shown by Fig. ^. Using the decompo- 

sition of the velocity field onto the different scales yielded 
by a MultiResolution Analysis with Daubechics' wavelets, 
using the scaling function (/)4( Daubcchies, 1992| ), we show 
that the correlation reaches ^ 0.9 at the largest scale avail- 
able. 

In Fig. ^, we plotted the actual velocity as a function 
of the measured velocity for various length scales. The 
clouds of points clearly show that granules motions sta- 
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Fig. 3. Dependence of the correlations with the spatial 
scale for the three time averages; the depth used is the 
optimal depth and the thickness is 78 km. The 2'^ scale 
corresponds to a resolution of 7x7 pixels or 667x667 km^; 
we understand the saturation of the 2^ scale as a conse- 
quence of the very few resolution elements left for this 
scale (~ 5x5). Line type have the same meaning as in 
Fig. I 



tistically underestimate the actual plasma velocity by a 
factor which (likely) tends to unity as the scale increases. 
When no filtering is made, measured velocities miss "real" 
velocities by roughly a factor 2. This quantitative dis- 
agreement is of course even more pronounced in derivative 
quantities (divergence or vorticity). 

This behaviour is of course no surprise because gran- 
ules are far from being passive lagrangian tracers: on the 
contrary they are active vortical structures which can 
move in the background fluid thanks to their own vor- 
ticity or the one of their neighbours. Their motion may be 
compared (but just qualitatively) to the random motion of 
molecules in a gas: only long time averages or large-scale 
averages are able to raise the signal of the mean motion 
above the noise of random motions. 

Finally, let us mention that we have done these tests 
using the two presently known methods of granule track- 
ing, namely LCT and CST. As illustrated in Fig. |l|, both 
methods give remarkably close results (within a few per- 
cent in correlation), a fact which gives confidence in the 
robustness of the results. CST, however, offers additional 
informations on the way individual granules follow the 
background flow; namely, we can appreciate which gran- 
ules are the most faithful tracers and characterize them by 
some property (size or lifetime for instance). We discuss 
this issue in the next section. 
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Fig. 5. Maximum correlation (solid line) between the ac- 
tual velocity field and the granule displacement at the 
place of the granule versus the surface area of the gran- 
ule. Large granules are only slightly more correlated than 
smaller ones. The dashed line indicates the depth at which 
correlation is optimum. 



3. Lagrange's view 

To see which granules are good or bad lagrangian tracers 
we computed the correlation between the mean velocity of 
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Fig. 6. Maximum correlation (solid line) between the ac- 
tual velocity field and the granule displacement at the 
place of the granule versus the lifetime of the granule. 
The dashed hnc indicates the depth at which correlation 
is optimum. 

individual granules (i.e. the velocities issued from granule 
tracking which yield the 45x45 dataset) and the actual ve- 
locity at the place of the granule. We did this computation 
for different granular sizes (Fig. ||) and lifetimes (Fig. 

Fig. H clearly shows that the size is a poor criterion 
for selecting granules whose motions represent the plasma 
flow. This figure, however, shows that large granules are 
sensitive to 'deep' undercurrents. The depth for optimal 
correlation increases with the size of the granules. 

On the other hand, the life-time is a good criterion for 
sampling the plasma velocity field. The motion of long- 
lived granules can reach 0.9 correlation with the actual 
flow field as shown in Fig. ^. Besides, the layer sampled by 
these granules is not precisely defined and oscillate around 
the T = 1 level. 

4. Conclusions 

We have used a simulation of compressible convection with 
radiative transfer in grey approximation to test the ability 
of granules at tracing the actual plasma flow. The box used 
for this simulation is 30x30 Mm^ wide and 3 Mm deep, 
resolved by a 315x315x82 grid. The results of these tests 
show that 

— Granules tend to be lagrangian tracers when the time 
and length scales of the flow tend to infinity: it shows 
that scale separation is a necessary condition for using 
granules at representing plasma flows. Quantitatively, 
we find that the length scale needs to be larger than 
2.5 Mm and the time scale longer than Ih for the cor- 
relation to be higher than 0.9. 

— They underestimate the velocity field, all the more that 
scale separation is weak. 



— Statistically, they probe a layer 300-400 km beneath 
the T = 1 surface. 

— Long-lived granules are good tracers. 

Hence, we see that granules can be used as tracers to 
reveal flows at meso- and supergranular scale. At smaller 
scales their own velocity field has a too strong interaction 
with the background velocity field. However, it may well 
be that velocity fields associated with ex ploding granules 



or wi th 'strong positive divergences' (see fiieutord et al 



2000 ), which are near the lowest (allowed) scale, are cor- 



rectly represented by granules motions, at least qualita- 
tively, since they are strong advective motions in nature. 

Finally, let us note that the situation in the real sun 
may not be better than that of the simulation for the 
Reynolds number is much higher: nonlinear interactions 
are indeed stronger and thus real granules are less pas- 
sive. On the other hand, the drift of the thermal structure 
with respect to pure advection, which is another pitfall 
of granule tracking, is likely correctly represented by the 
simulation since it depends on the Peclet number which 
is modeled accurately by the simulation in the layers of 
interest. 
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